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Abstract: Due to the strong and steady TeV 7-ray emission from the Crab Nebula supernova remnant, its 
measured flux and energy spectrum can be used to verify the calibration and data reduction methods ap- 
plied to IACT data acquired over many observing seasons. This gives us confidence in the results obtained 
on variable TeV sources observed over the same period and in relating the sensitivity of new instruments 
to historical datasets. Here we present the results of an analysis of 65.3 hours of good quality data taken 
on the Crab Nebula between October 2000 and March 2006 with the Whipple 10m telescope. The total 
exposure resulted in a 46 a signal with 1 1886 selected excess events. The energy spectrum was best fit by 
a power law of the form dN/dE = (3. 19 ± 0.07 stat .) x 10~ n ■ (E/lTeV)" 2ti4±0 03stat cm" 2 s _1 TeV _1 
in the energy range 0.49-8 TeV. The systematic uncertainty in the flux was estimated to be 30%, with a 
systematic error of 0.2 in the photon index. A reasonable agreement is shown for a fit to a constant flux 
over the 6 years. 



Introduction 

The Crab Nebula supernova remnant has served as 
the TeV 7-ray "standard candle" for Imaging At- 
mospheric Cherenkov Telescopes (IACTs) since its 
successful detection with a 37 pixel camera on the 
Whipple 10m telescope in 1989 [1,2,3,4]. Here 
we provide a record of its flux and energy spectrum 
from 2000 to 2006 as observed with the Whipple 
10m telescope in its current configuration. 

The aim of this work is to prove the suitability of 
various data reduction and analysis techniques and 
the stability of the instrument for long-term source 
monitoring. The full study will appear in [5]. Here 
we compare the results obtained with two sets of 
Hillas image parameter-based background rejec- 
tion cuts, "Hard cuts" and "Loose cuts". We dis- 
cuss errors in energy reconstruction and the varia- 
tion in effective collection area with the zenith an- 
gle of observations, both of which were quantified 
using simulations. Finally, we plot three integral 
flux measurements per observing season as a Crab 
Nebula light curve, and a differential energy spec- 
trum obtained from the combined data set is pre- 



sented alongside results from HEGRA, H.E.S.S. 
and MAGIC. 



Observations 



Obs. Period 


T L 


Nruiis 


(©> 


(Rate) 


(Month/Year) 


(hr) 




(cleg.) 


(Hz) 


10/00 - 01/01 


13.2 


29 


18.8 


29.7 


10/01-03/02 


23.9 


52 


17.1 


28.0 


11/02-02/03 


8.3 


18 


17.2 


17.6 


11/03-03/04 


5.5 


12 


17.8 


23.5 


10/04-04/05 


5.5 


12 


15.0 


22.0 


10/05 - 03/06 


8.8 


19 


18.6 


22.8 


Total 


65.3 


142 


17.6 


25.5 



Table 1: Observations of the Crab Nebula. Tl is 
the total time ON source, (Rate) is the mean raw 
trigger rate. 



The Whipple 10 m telescope is operating with a 
379 pixel, 2.6° field of view camera. This is the 
high-resolution inner section of the 490 pixel cam- 
era installed in 1999 [6]; the 111 outer guard ring 



Crab Nebula with the Whipple 10 m 
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Figure 1: Mean relative error in reconstructed 
7-ray energy as a function of simulated energy 
Ettuo at 20°, 30° ,and 40° zenith angle. 



pixels were decommissioned in 2003 and conse- 
quently have not been included in this analysis. 

The Crab Nebula datasets used are listed by sea- 
son in Table 1. For this systematic study only 
28 minute observation runs centered on the source 
position (ON runs) followed by a matching OFF 
source run were included. The data cover a range 
in zenith angle of 10-30°. Only data taken un- 
der good weather conditions, for which the RMS 
spread in the raw telescope trigger rate was less 
than 1.5 Hz were included. As can be seen in Ta- 
ble 1, the average raw trigger rate varies somewhat 
between observing seasons, which is indicative of 
small changes in telescope efficiency probably due 
to a combination of instrumental and environmen- 
tal factors. 

Data Reduction 

We applied the Islands method of [7] to extract 
clean Cherenkov images. The images were param- 
eterized according to image intensity, shape and 
orientation. The Length and Width were converted 
to the "reduced scaled" parameters RSL and RSW 
[2]- 

A subset of 10 hours of the Crab Nebula obser- 
vations from 2000-2006 recorded at ~20° zenith 
angle were used to select two sets of 7-ray selec- 
tion cuts: Loose cuts and Hard cuts, the latter for 
measuring a 7-ray flux and energy spectrum with 
a high detection significance (following equation 
17 of [8]). The number of excess 7-ray type events 
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Figure 2: Effective collection area versus true sim- 
ulated energy after Hard cuts. 



N 7 = Nqn — a • Nqff is calculated from number 
of events passing the cuts, scaled by the ratio of the 
ON to OFF source exposure time a. Table 2 lists 
the multi-parameter selection cuts, corresponding 
significance and 7-ray detection rate (defined as 
R = N 7 / Ton). 



Energy Evaluation and Effective Area 

Simulated 7-rays at four zenith angles and three 
telescope efficiencies /1 were used to fill lookup ta- 
bles for the mean simulated energy Etiuc, Length, 
and Width as a function of Distance and Log(Size). 
A two-dimensional Gaussian smoothing function 
was applied to the lookup tables with o"Dist = 0.05° 
in Distance and 0Log(s) = 0-01 in Log(Size). For 
each event, the reconstructed energy Er oc was cal- 
culated from the lookup tables by linearly inter- 
polating between /1 and cos <d. Figure 1 shows 
the mean relative error of the reconstructed energy 
(ER ec — E-prue) / Eiruc as a function of simulated 
energy E<i> ue at zenith angles of 20°, 30°, and 40° 
after applying Hard selection cuts. At low ener- 
gies, ER ec is overestimated due to events with in- 
tensity near the telescope trigger threshold. A us- 
able energy range above Es a fc with a relative error 
of < 10% was determined for each zenith angle, 
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Cut 


Size 
(pe) 


Distance RSL 
(deg.) 


RSW 


Alpha 
(deg.) 


Len./Size 
(deg./pe) 


cr/Vhr 


R-y 

(min. -1 ) 


Hard 

Loose 


>80 
>80 


0.2-0.95 -2.0-1.6 
0.2-0.95 -2.0-2.0 


-2.0-1.6 
-2.0-2.0 


<15 
<22 


<0.0011 


5.26 
3.56 


3.12 ± 0.18 
5.51 ±0.49 



Table 2: Image selection cuts. The significance a/Vhr and rate R 7 of selected events are from 10 hours of 
observations of the Crab Nebula during 2000-2006 at ^20° zenith angle. 
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Figure 3: Contour plot of the 68%, 95% and 99.9% 
confidence intervals from the \ 2 fit to a power law 
for the total 2000-2006 dataset after applying Hard 
cuts or Loose cuts. 



with the minimum being Es a fc = 0.5 TeV at zenith 
angle 20°. 

In order to account for biases in energy reconstruc- 
tion, the energy spectrum is measured using the ef- 
fective area as a function of reconstructed energy. 
The maximum distance R = 400 m of the sim- 
ulated air showers from the telescope was chosen 
to encompass the full impact parameter range of 
triggered events (0-270 m). The total number of 
simulated 7-rays is represented by N s i m .(E,0,/i), 
and the number of detected events passing selec- 
tion cuts as N so i.(E,0,/i). The effective areas were 
fitted with an analytical function modified from 
equation 3 of [9], 

Figure 2 shows the effective areas after Hard cuts 
as a function of true simulated energy Etiuc at 
zenith angles 6 of 20°, 30°, and 40° with telescope 
efficiency 85%. 
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Figure 4: Total 2000-2006 measurement of the 
Crab Nebula energy spectrum compared to previ- 
ous measurements. 



The integral 7-ray flux above a chosen energy 
threshold is: 



'>E th 



/ dF \ 

(i-r) 



E 



(2-r) 

th 



where a fixed source spectrum with power law pho- 
ton index T is assumed. The excess number of 
events N^ 

N, 



-. 7 is given by 

dF r r 
dE Je 



/„ L Aff(E,e(t),/i(t))-E- r dtdE 

The normalisation factor (^§)th can thus be esti- 
mated from the N 7 measured for each run and the 
integral of the differential rate multiplied by the to- 
tal livetime of the run Tl- 



Results and Conclusions 

A total Crab Nebula energy spectrum over 2000- 
2006 was measured with both Hard and Loose se- 
lection cuts. Figure 3 shows a contour plot from 
the x 2 fit errors in flux normalization factor ( j§)th 



Crab Nebula with the Whipple 10 m 



Dataset 


a/^/hr 


KdEh TeV 


r 


(Year) 








'00- '01 


5.32 


3.38 ±0.19 


2.57 ±0.10 


'01 - '02 


6.19 


3.10 =h 0.11 


2.54 ± 0.06 


'02 - '03 


6.91 


2.75 ±0.16 


2.66 ±0.11 


'03 - '04 


5.21 


3.75 ± 0.29 


2.45 ±0.10 


'04 - '05 


5.54 


2.90 ± 0.26 


2.68 ±0.19 


'05 - '06 


5.28 


3.58 ± 0.25 


2.60 ±0.12 


Tot-Hard 


5.64 


3.19 ±0.07 


2.64 ± 0.03 


Tot. Loose 


3.48 


3.31 ±0.10 


2.62 ± 0.04 



=~2000 ; 2001 



Table 3: Energy spectrum of the Crab Nebula with 
Hard cuts. The power law flux normalization 
at 1 TeV is in units of 10~ n cm" 2 s" 1 TeV" 1 . 



and photon index T for the total Whipple 10 m 
2000-2006 dataset with Hard and Loose cuts. The 
best fit values agree at the 2a level. Using Hard 
cuts, the best fit model was a power law over the 
energy range 0.49-8 TeV with: 



dN 
dE 



(3.19 ±0.07) x 10~ n ■ ( 



E 



ITeV 



_1 TeV" 



Figure 4 shows the corresponding total energy 
spectrum of the Crab Nebula compared to past 
measurements with the Whipple 10 m in 1994- 
1995, HEGRA in 1997-2002, and H.E.S.S. in 
2003-2005 [10] [2] [4]. 

Table 3 lists the results from a power law fit to the 
energy spectrum for the individual observing sea- 
sons using the Hard cuts (a more comprehensive 
table including Loose cuts results can be found in 
[5]). A reasonable agreement in the measured pho- 
ton index T and flux normalization factor (4fr)th 
is found between data sets within statistical er- 
rors. The RMS spread in T between datasets is 
0.06, with a mean statistical error of 0. 1 1 from each 
dataset. 

The integral flux F>it c v was calculated from the 
fitted power law spectrum for each dataset. As can 
be seen in figure 5, the measured values were con- 
sistent with steady emission over 6 years; between 
datasets the RMS spread in the integral flux was 
12%. This gives us a sound basis for the character- 
isation of the flux variability of other sources. 
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Figure 5: Integral flux of the Crab Nebula over 
time. 
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